Stability of Fe-based alloys with structure type CeCr 2 3 
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Bulk metallic glass forms when liquid metal alloys solidify without crystalization. In the search 
for Iron-based bulk glass- forming alloys of the metal-metalloid type (Fe-B- and Fe-C-based) , crystals 
based on the structural prototype CeCr23 often preempt the amorphous phase. Destabilizing this 
competing crystal structure could enhance glass-formability. We carry out first-principles total 
energy calculations of enthalpy of formation to identify third elements that can effectively destabilize 
CeCr23. Yttrium appears optimal among transition metals, and rare earths also are suitable. Atomic 
size is the dominant factor. 

PACS numbers: 61.50.Lt,61.43.Dq, 71.20.Be, 81.30.Bx 
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I. INTRODUCTION 

Iron-based amorphous alloys are used in transformer 
cores, where their low magnetic coercivity reduces en- 
ergy loss. Popular glass-forming alloys are based on Fe 
together with metalloid elements such as B or C. Bulk 
Iron-based amorphous alloys could become important 
structural materials, but optimal glass-forming compo- 
sitions are not yet known. Multicomponent alloys con- 
taining fourth row transition metals and rare earths show 
promise 0,Q,I3- 

We previously explored the quaternary B-Fe-Y-Zr 
phase diagram [j] , identifying stable and metastable crys- 
tal phases and computing their enthalpies of formation. 
This study identified crystalline structures based on the 
CeCi'23 prototype as important competitors to glass for- 
mation. It appeared that the competition is more prob- 
lematic in the case of B-Fe-Zr than in the case of B-Fe-Y. 
To ensure the optimal selection of alloy system, we now 
carry out a systematic study of many candidate "third el- 
ements" and compare them with regard to stability of the 
CeCr23 structure. We do this both for the case of B-Fe- 
and C-Fe-based alloys. We show that atomic size mis- 
match destabilizes the C 6 Cr 2 3 structure for sufficiently 
large atoms such as Yttrium and rare earths [EE3] 



II. METHODS 

Our ab-initio calculations use the program VASP 
together with the projector-augmented wave method, 
an all-electron generalization of the pseudopotential ap- 
proach 0, ^} . We employ the Perdew-Wang generalized 
gradient approximation 11] (GGA) exchange-correlation 
functional with the Vosko- Wilkes Nussair jlj] spin inter- 
polation. These choices give excellent results for bulk el- 
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emental Fc 10]. GGA is needed instead of LDA is neces- 
sary to properly reproduce magnetization and lattice con- 
stants |13| ■ Our magnetic calculations are spin-polarized 
(i.e. collinear magnetization) and are employed for any 
structure containing 50% Fe or higher. All calculations 
for Carbon-based binaries and ternaries are performed 
at a constant cutoff energy of 400 eV, the default for our 
Carbon potential. All calculations for Boron-based bina- 
ries and ternaries are performed at 320 eV, the default 
for our Boron potential. More details and discussion of 
convergence, etc. are given in Ref. Q. All the ab-initio 
data on which this paper is based can be obtained on the 
WWW at Ref. 0. 

The composition space of an A^-component alloy is a 
set of N composition variables forming aniV-1 dimen- 
sional simplex (respectively, a point, line segment, trian- 
gle and tetrahedron for N — 1,2,3,4). Structural en- 
ergies form a scatter-plot over this simplex. Stable low 
temperature phases lie on vertices of this scatter-plot. 
Edges and facets of the convex hull represent coexistence 
regions of the phases at adjoining vertices. 

The tie-lines and tie-planes connecting all pure ele- 
ments in their lowest energy structures forms a useful 
reference for alloy energies. The distance AHf or of an 
alloy energy from the tie-surface joining pure elements 
is known as its enthalpy of formation. It is an enthalpy 
because volume relaxation means we work at fixed pres- 
sure, P = 0. Strong compound formation is reflected in 
large negative enthalpy of formation. The value of AH is 
determined solely by the cohesive energy of a given struc- 
ture relative to the cohesive energies of its constituent 
pure elements. 

High temperature phases should lie above the con- 
vex hull, but be sufficiently close that entropic effects 
(e.g. atomic vibrations, vacancies or chemical substitu- 
tion) can stabilize them. Metastable phases also should 
lie close to the convex hull, so that their free energy is 
less than the liquid free energy at temperatures below 
freezing. Although AHf or is usually negative for high 
temperature and metastable phases, their energy differ- 
ence AE from the convex hull is small and positive. AE 
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FIG. 1: Enthalpies of formation and their convex hulls for the B-Fe (top) and C-Fe (bottom) binary alloys. Notation: heavy 
circles denote known low temperature phases, light circles denote known high temperature phases, diamonds denote known 
metastable phases, squares denote unreported structures. 



measures the thermodynamic driving force for decompo- 
sition into the appropriate combination of stable phases. 
In contrast to AH, the value of AE depends on the co- 
hesive energies of other competing structures. Discovery 
of a new stable structure will increase the assessed AE 
values of previously known structures. 

In principle all possible combinations and arrange- 
ments of atoms should be considered to ensure the op- 
timal possibilities are found. This is clearly impossible. 
Rather, we choose plausible structures for consideration 
by chemical substitution into known structures of simi- 
lar compounds [15J . We especially consider alloy systems 



with similar atomic size ratios or other chemical proper- 
ties. 

Structures are denoted using their prototype names 
and Pearson symbols. For example, we will be inter- 
ested in the CeCr23 prototype. The element Cr will be 
replaced by Fe, and in some cases the element C will be 
replaced by B. The Pearson symbol for the CgCr 2 3 proto- 
type is cF116, indicating cubic symmetry, face-centering, 
with 116 atoms per unit cell. The primitive cell of 
C 6 Cr 23 .cF116 contains 116/4=29=23+6 atoms. 

Using these methods, we built a database of structural 
energies. For a given .ZV-component alloy system of inter- 
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FIG. 2: Enthalpies of formation and their convex hull for the C-Y binary alloys. Notation as in Fig. 



est we extract from our database energies of structures 
containing all, of just some, of the chosen elements. We 
use a standard convex hull program (qhull 0) to iden- 
tify stable structures and the coexistence regions that 
connect them. Based on the output of this program, we 
calculate values of AHf or and AE for every structure. 
Numerical data for the compounds considered here, and 
more than 1500 others, can be found at Ref. 1141 . 



III. RESULTS 
A. Binaries 

Cohesive energies of binary B-Fe and C-Fe alloys are 
shown in Fig. ^ Out results for both alloy systems are 
in perfect qualitative agreement with experiment |!7lll8| . 
because all the known low temperature stable phases oc- 
cur on the convex hulls (notice C-Fe has no stable com- 
pounds) and all the know high temperature, metastable 
and unknown hypothetical structures lie above the hull. 
It is impressive how sensitive the density functional the- 
ory is to differences in chemical identity, with the many 
distinctions between B-Fe and C-Fe all being faithfully 
reproduced. 

Many metastable phases are known in the C-Fe binary 
system. Two of these, CFe4.cP5 and CFe3.hP8 are based 
on FCC Iron with Carbon intcrstitials in, respectively, 
tetrahedral and octahedral sites. The CFe3.oP16 struc- 
ture is also an important metstable phase in B-Fe. It can 
be considered as a strong distortion of the Fe3Si.cF16 
structure, caused by the small atomic sizes of B and C 
relative to Si. Similarly, the energies of the interstitial 
cP5 and hP8 structures are lower for C than for B be- 



cause the Carbon atoms are smaller than Boron, 
chief B-Fe structures are discussed further in Ref. 



The 



The lowest-lying metastable C-Fe binary compound 
has the structure of CgC^. This crystal structure ap- 
pears to be the most important competitor to metallic 
glass formation, and further destabilizing it is the goal of 
this work. 

To evaluate stabilities of a ternary alloy system such 
as C-Fe-Y we need to examine all three of its binary 
subsystems. The Fe-Y diagram was previously discussed 
in |j] and need not be repeated here. The C-Y binary 
phase diagram is poorly known |l9ll2fjj |. with three phases 
(/3-CigYi5 and a, /3-C3Y2) of unknown structure. The 
reported stable low temperature tP68 structure of a- 
C19Y15 (isostructural with CigSci5.tP68) suffers from 
large displacements during relaxation, and a final re- 
laxed energ y ab ove the convex hull. We find instead that 
C4Y3.tP70 [I3 is stable at low temperatures. 

The high temperature 7 phase probably takes the re- 
ported Fe4N.cF8 structure in the Y-rich limit. This 
structure is based upon an FCC lattice of Y atoms with 
C occupying octahedral interstitial sites. When these 
sites are fully occupied the unit cell composition is C4Y4 
and the structure becomes NaCl. On the C-rich side, 
we begin occupying tetrahedral interstitials and the en- 
ergy rapidly grows beyond values that are plausible for a 
high temperature phase. We also tested Y vacancies and 
found them even higher in energy. We believe the 7 phase 
actually terminates at 50% Carbon, rather than the re- 
ported 67%. The postulated /3-C 2 Y.cF12 structure 
at 67% Carbon (based on an FCC lattice of Y with C 
fully occupying all tetrahedral interstitials) is highly un- 
stable and presumably incorrect. 
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FIG. 3: Energy diagram for ternary C-Fe-Y. Plotting symbols as in Fig.Q Values in parenthesis indicate enthalpy of formation 
AH f or for stable compounds, and energies above the convex hull AE for metastable structures. 



B. Ternaries 



Having verified the ability to reproduce these binary 
phase diagrams, we now investigate ternary systems. 
Fig |3 illustrates a typical ternary energy diagram, in 
this case for the alloy C-Fe-Y. This diagram confirms 
stability of previously known ternary phases and also 
proposes B2Fe2Y.tI10 as stable, though this structure 
has not previously been reported |21| . We also propose 
C2FeSc.tP8 as the probable structure of a previously re- 
ported metastable phase of unknown structure. Similar 
energy diagrams for all ternaries discussed below may be 
found at Ref. 0. 

We now focus our attention on one specific structure, 
CeCr23, and examine its stability for a variety of B-Fe- 



and C-Fe-based ternaries. Although Chromium-based in 
the prototype structure, it is a well known metastable 
phase in Iron-based alloys such as C-Fe-Mo and we pre- 
dict it is actually stable in B-Fe-Sc, B-Fe-Nb and B-Fe- 
Mo. 

Our plan is to find elements that can mix with the bi- 
naries in the liquid state but that will destabilize CgCr23 
due to their large atomic size. The requirement of misci- 
bility in the liquid state rules out alkali metals and alkali 
earths. The requirement of large atomic size rules out 
middle and late transition metals. We thus identify the 
primary candidates as the early transition metals and the 
rare earths. 

Our discussion begins with Yttrium, the first transition 
metal of row 4 in the periodic table. The first question to 
address is the optimal site for Y atoms. As a large atom, 
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TABLE I: Site energies (AE (eV/atom)), near-neighbor dis- 
tances and Voronoi volumes for Y substitution on different 
Wyckoff classes. 



site 


4a 


8c 


32f 


48h 


energy 

TV 

V 


0.135 
2.62 
20.1 


0.085 
2.71 
24.5 


0.145 
2.42 
23.7 


0.146 
2.40 
22.4 



it cannot enter as an interstitial, nor can it substitute for 
Carbon. There are thus 4 plausible sites, the Iron sites 
of Wyckoff classes 4a, 8c, 32f or 48h (see Table 0. 

Table H] evaluates the relaxed energies of a single 29 
atom primitive cell of cF116 with one Y atom replacing 
one Fe atom, respectively, on each distinct Fe site. Ev- 
idently Y atoms favor Wyckoff class 8c, most likely as 
a result of their large size. Indeed, placing Y on site 8c 
results in both the largest Voronoi volume V for Y atoms 
and the largest near- neighbor distance ry- 

We also checked that the 8c site is prefcrcd for a sec- 
ond Y substitution, after the first Y atom is already on 
8c. Since there are only two 8c sites in a primitive cell 
(8/4=2) this corresponds to complete filling of the 8c 
sites with Y atoms. Double occupancy of the 8c site 
by the large atom is consistent with the solved crystal- 
lographic structure of BgCo2iZr2, which is isostructural 
with CeCr23. 

The following energy study focuses on the case of dou- 
ble substitution to fully occupy the 8c sites. The en- 
ergy cost for occupying both 8c sites by Y is more than 
twice the cost for occupying just one, as can be seen 
by comparing tables [I] and [H] The energies for 0, 1 
and 2 Y-atom substitutions on site 8c in C6Fe23 are, 
respectively, A£=0.043, 0.085 and 0.141 eV/atom in a 
29 atom primitive cell. So the first substitution costs 
29 x (0.085 - 0.043) = 1.22 eV, while the second costs 
29 x (0.141 - 0.085) = 1.62 eV. Presumably this is be- 
cause lattice strains created by insertion of the first Y 
atom are accomodated partly by shrinkage of the volume 
around the second 8c position. When the second Y atom 
is introduced this accomodation is no longer possible. 

A picture of the cF116 structure at composition 
CeFe2iY2 is shown in Fig. The 8c sites have coor- 
dination number Z = 16, with only Fe atom neighbors. 
The Voronoi polyhedron (Watson and Bennett |22| class 
(0,0,12,4)) has 12 pentagonal faces and 4 hexagonal faces, 
characteristic of local tetrahedral close packing. The 
Voronoi polyhedron is illustrated in Fig.0] However, the 
remainder of the structure does not exhibit tetrahedral 
packing. 

The Carbon sites, Wyckoff class 24e, have a Voronoi 
polyhedron (Watson and Bennett class (0,8,0)) with 8 
square faces and a highly nongeneric four-fold vertex. 
The eight neighbors of each Carbon atom are Iron atoms 
at distances of 2.05 and 2.15 A, substantially greater than 
the 1.79-1.99 A found in other metastable C-Fe struc- 
tures. This loose binding of the Carbon atom probably 
explains why the CgC^ structure is actually more fa- 



TABLE II: Energies for complete TM substitution on site 8c. 
the top row are B-Fe and C-Fe binaries for comparison. 



z Ch 


B-Fe 


C-Fe 




A TP 


A TJ 
l\tl 


rr 


A TP 


A TJ 
l\tl 


Tt 


re 


0.018 


-0.168 


2.40 


0.043 


0.043 


2.42 


39 y 


0.074 


-0.163 


2.61 


0.141 


0.059 


2.61 


40 Zr 


0.016 


-0.253 


2.54 


0.025 


-0.037 


2.54 


41 Nb 


stable 


-0.246 


2.48 


0.002 


-0.028 


2.48 


42 Mo 


stable 


-0.212 


2.42 


0.030 


0.022 


2.41 


' 2i Sc 


stable 


-0.244 


2.53 


0.054 


-0.035 


2.53 


57 La 


0.094 


-0.092 


2.60 


0.195 


0.189 


2.61 


89 Ac 


0.232 


0.046 


2.66 


0.228 


0.228 


2.65 



vorable for B-Fe than for C-Fe, because Boron is slightly 
larger. 

Now we compare the stabilities with different choices of 
large atom. Our data is presented in Table ITT1 for several 
choices of large atom. In each case we verified that the 
8c site is prefered. We fully occupy the 8c site, with 
two large atoms per face-centered primitive cell. Values 
of tt are distances from large atom to nearest neighbor 
transition metal in the fully relaxed structure. 

Evidently, of the 4d transition metals, the largest 
atom, Yttrium, is the most effective at destabilizing 
CeCr23. Indeed, Zr, Nb and Mo tend to stabilize it rela- 
tive to the B-Fe and C-Fe binaries. However, too large a 
Y concentration is dangerous for glass formation, because 
of the low-energy metastable C3Fei7Y2.hR22 structure 
(Y content 9.1%) and the stable structure CFei 4 Y 2 .tP68 
(Y content 11.8%) visible in Fig.EJ 

We can also test whether the 4d row is an optimal 
choice by examining other members of group IIIA (the 
first column of transition metals, consisting of Sc, Y, La 
and Ac). Clearly the 3d element Sc is too small to desta- 
bilize CgCr23. Lanthanum (5d) and Actinium (6d) are 
both highly effective. The lanthanide and actinide rare 
earth elements, which behave chemically like Lanthanum 
and Actinium, should likewise do well. Of course, Ac- 
tinides may present other difficulties not considered here! 



IV. CONCLUSIONS 

We compared the effectiveness of different large-atom 
substitution as a way to destabilize the CgG^ struc- 
ture of B-Fe and C-Fe. The first task was to identify the 
most energetically favorable substitution site for the large 
atom. We found, in every case considered that Wyck- 
off site 8c (a locally tetrahedrally close-packed site) was 
most favorable. We showed that among the 4d transition 
metals Yttrium was the only element capable of destabi- 
lizing the CgCr23 structure relative to the binary alloy. 
Comparing group IIIA elements, we found that Sc was 
inadequate but Y, La and Ac are fine. 

Clearly atomic size is a crucial consideration for desta- 
bilizing this structure. Indeed, the 8c site prefers atoms 
slightly larger than Fe, so replacement with Mo, Nb or 
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FIG. 4: (left) Structure of C6Fe21Y2.cF116 viewed along (1,1,1). Atom size indicates height. Hash marks are 1 A. (right) 
Voronoi polyhedron of Y site, Wyckoff class 8c, viewed along (1,0,0). 



Sc can actually stabilize cF116. We conclude that Yt- 
trium and rare earth elements can enhance glass forma- 
tion. However, too large an Yttrium concentration (9% 
or more) can lead to formation of other competing crystal 
phases such as C 3 Fei 7 Y 2 .hR22 or BFei 4 Nd 2 .tP68. Al- 
kali earth elements such as Ca are also highly effective 
at destabilizing CgC^, but they do not mix well with 
Fe-based alloys, at least at atmospheric pressure. 
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